Glassy Motion of Elastic Manifolds 
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We discuss the low-temperature dynamics of an elastic manifold driven through a random medium. 
For driving forces well below the T = depinning force, the medium advances via thermally activated 
hops over the energy barriers separating favorable metastable states. We show that the distribution 
of waiting times for these hopping processes scales as a power-law. This power-law distribution 
naturally yields a nonlinear glassy response for the driven medium, v ~ exp(— const x 



On 
ON 

u 

Or 

< 

O 

CO 



> 

oo 

o 



c 

o 
o 



The dynamics of elastic manifolds in random media has 
received an intense attention during the last decade and 
is in a state of rapid current development. The reasons 
for this interest are twofold. On one hand, elastic man- 
ifolds in a random environment are one of the simplest 
examples of a glassy system and yet exhibit a very rich 
static and dynamical behavior. Secondly, they can be 
used as models to study both nonlinear collective trans- 
port in driven disordered systems, such as charge density 
waves, polymers, driven interfaces, dislocations in solids, 
and magnetic flux lines in type II superconductors, as 
well as stochastic kinetic processes, such as stochastic 
growth and kinetic roughening (see (TJQ] for a review). 

The dynamics of driven elastic manifolds is the result 
of the interplay between quenched disorder and the inter- 
action among the many elastic degrees of freedom that 
compose the manifold. A key physical quantity is the av- 
erage velocity v of the driven manifold as a function of the 
applied force F. At zero temperature there is a depinning 
transition from a pinned state where v — to a sliding 
state at a critical driving force F c . A finite temperatures 
washes out the transition and the mean velocity is then 
finite for all driving forces. For low temperatures and 
driving forces well below the T = depinning threshold, 
F c , the dynamics of the driven manifold is controlled by 
thermally activated jumps of correlated regions over the 
pinning energy barriers separating different metastable 
states. In this region the mean velocity is highly non- 
linear and has been evaluated via a scaling approach , 
with the result, 



v~exp[-U(F)/T}, 



(1) 



where U(F) is the optimal energy barrier for creep be- 
tween favorable metastable states. Under the action of 
the external drive, sections of the manifold that are ini- 
tially pinned move to a more favorable metastable state, 
determined by the condition that the energy gain due 
to the driving force equals the elastic deformation and 
pinning energy of the medium. For F << F c large sec- 
tions of the manifold hop long distances to find the next 
optimal-energy state. This yields a large energy barrier 



that diverges algebraically with vanishing driving force, 
U(F) ~ (F c /F)^, with fi a characteristic exponent. The 
form given in Eq. |l| was first proposed for the motion 
of dislocations in crystals || and is now widely used to 
describe the low temperature dynamics of vortex lines in 
type-II superconductors pj. 

An important recent development has been the confir- 
mation by numerical studies of an elastic string in d = 2 
of the basic assumption of the scaling approach that 
the barriers between the metastable states scale with the 
length of the fluctuating string segment in the same way 
as the fluctuations in the free energy (aside from loga- 
rithmic corrections) 0. Despite a substantial analyti- 
cal Jl-Q] and numerical |§-|l3|] study of the dynamics of 
driven disordered elastic manifolds, many open questions 
remain. 

In this Letter we focus on the nature of the mecha- 
nism by which the manifold selects the appropriate op- 
timal segment controlling the dynamics. For simplic- 
ity we address this question for the case of an elastic 
string (magnetic flux line or interface) driven through 
two-dimensional point disorder. The paper is organized 
as follows. First we briefly review the results of the scal- 
ing theory of the thermally activated motion of elastic 
manifolds in disordered media following || . We then de- 
velop a description of the avalanche-like string motion. 

The Hamiltonian of a D-dimensional elastic manifold 
driven through a d-dimensional disordered medium is 



H = J d d x 



C fdu\ 



u 



(2) 



where C denotes an elastic stiffness constant (e.g., the 
linear tension of the string; for simplicity we consider 
the isotropic medium) and u(x, t) is the n-dimensional 
transverse displacement field of the manifold. Eq. (|^) 
describes for instance an elastic string (D = 1) in two 
(d = 2, n = 1) and three [d = 3, n = 2) dimensions, and 
two- or three-dimensional vortex lattices (d = 2, D = 2, 
n = 2 and d = 3, D = 3, n = 2, respectively). The 
disorder consists of uncorrelated point defects of density 
m. It is described by a random potential U (x, u) of 
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range £ and variance A = voy/rii^ d , with vq the maxi- 
mum depth of the potential well of a single pinning cen- 
ter. In the absence of driving force the string adjusts to 
the random landscape and traverses the medium along 
rough optimal paths determined by balancing the elas- 
tic and pinning energies. The geometry of these opti- 
mal paths is characterized by the roughness of the man- 
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ifold, defined as w(L) = ([u(x + L) — u(x)] ) , where 
(...) denotes the average over both thermal fluctuations 
and quenched disorder. At large distances the roughness 
scales as(i)~( (L/L c ) , where C < 1 is the roughness ex- 
ponent and L c is the pinning correlation length. We con- 
sider the case of weak disorder such that L c 3> £. In this 

regime domains of linear size L c = £ (C£ D /A) 2 °^ are 
pinned coherently by disorder when F — 0. The pinning 
length L c is the smallest scale on which barriers between 
metastable states exist at F = 0. The minimum average 
energy barrier between neighboring metastable positions 
of a pinned segment L c is U c — C^ 2 /L C D . Optimal 
metastable configurations of sizes L > L c are then sep- 
arated by barriers U(L) ~ U C (L/L C ) 2<:+D - 2 . At T = 
the string starts to slide when the applied force F can 
depin a region of linear size L C1 yielding a threshold force 
for sliding F c = C£/L 2 C . 

We now consider the dynamics at a finite but low tem- 
perature T <C U c , where the elementary pinning scales 
U c and L c are not renormalized significantly by thermal 
fluctuations (lj. Under the action of the driving force a 
domain of size L can be displaced to a new more favor- 
able metastable state. For F <C F c the dynamics can be 
described as the nucleation of an elementary excitation 
or nucleus of size L. The free energy cost for creating 
such a nucleus is given by jF[L] = U C {L/ L C ) 2C,+D ~ 2 — 
FL c J £ i (L/L c )'' +D . The size L opt of the optimal excita- 
tion is obtained by minimizing this free energy cost, with 
the result L opt = L c (F c / F) 1/{2 ~° . Nuclei with L < L opt 
collapse, while nuclei of size L > L opt expand. In other 
words excitations on scales larger than L opt slide freely, 
while pinning barriers on scales L < L opt can be over- 
come only via thermally activated hops. The length scale 
L opt determines the upper bound above which thermally 
activated processes are no longer relevant. The creep rate 
is determined by the corresponding optimal energy bar- 
rier, U(F) = T(L opt ), with U(F) = U C {F C /Ff, where 
/i = (2£ + D — 2)/(2 — (). For an elastic string in two 
dimensions £ = 2/3, D = 1, and fi=l/4. 

We now show that the form given in Eq. ([[]) for the 
mean motion arises naturally when the driven manifold 
sequentially nucleates all possible excitations on length 
scales L < L opt to overcome all the possible energy bar- 
riers separating metastable states on scales E < U(F). 
For the sake of simplicity we consider an elastic string 
in 2d. We demonstrate that the distribution of waiting 
times t(E) for hops between metastable states separated 
by energy barriers E < U (F) scales as a power-law, 



*(t)~t-'-<*, (3) 

with a < 1. This microscopic mechanism of string dy- 
namics can be viewed as an avalanche-like motion. The 
distribution of waiting times given in (g) is cutoff at the 
waiting time t[U(F)] corresponding to the optimal bar- 
rier U(F). For hopping times r exceeding t[U(F)] only 
very rare barriers with E 3> U(F) can retard the motion 
(rare events) and the distribution function decays much 
faster than in Eq. (g). The average waiting time that 
controls the mean string velocity is therefore determined 
by the cutoff r[U(F)]. 

We begin by considering a segment of finite length 
L < L opt initially pinned in a metastable configuration. 
We derive an expression for the probability Pl{E) that 
this segment encounters energy barriers smaller than E 
when it samples all possible departures from its initial 
configuration under the action of the external force F. 
For any finite driving force the set of all configurations 
that the segment may sample during its activated motion 
form a connected cluster of temporarily pinned states. To 
describe this cluster, we fix the ends of the segment L 
and imagine dividing it in elementary units of longitu- 
dinal size L c . An elementary move consists of the hop 
of the i-th unit L c across a transverse distance £ to a 
new metastable position by overcoming the elementary 
pinning energy barrier Ui. Here the {£/j}'s are random 
variables fluctuating around their average value, the min- 
imal average pinning barrier, U c . The segment L can 
be thought of as a singly-connected necklace of n (unit) 
beads, where the beads corresponds to the elementary 
sections L c . To begin, we consider the simplest class of 
departures of L from its initial configuration and describe 
the advance of L via the hopping of individual beads over 
the fluctuating barriers Ui. The barrier associated with 
the motion of the segment L is then the largest of the 
barriers of the elementary hops L c . To understand this 
we note that the hop of a section L of the string to a 
nearby metastable state can be thought of as the motion 
of the domain boundary between spin up and spin down 
regions in an Ising spin system with random couplings 
{Ji}. Each elementary unit L c corresponds to a single 
Ising spin. The hop of the i-th segment L c over the cor- 
responding energy barrier Ui corresponds to a single spin 
flip. The flip of spin io along the domain boundary will 
cost an energy of order Ji ~i + Ji +i, where Ji -i and 
Ji a +i are the couplings to the neighboring spins. The 
flipping of this particular spin will facilitate the flipping 
of the neighboring ones, since now overturning spin io — 1 
will only cost an energy Ji -2, which may be larger or 
smaller than Jj _i + Ji +i- In either case the process 
of overturning the entire domain boundary will be con- 
trolled by the largest energy cost for flipping a single 
spin. If L is composed of n elementary segments L c with 
energy barriers Ui randomly fluctuating about U c , then 
the energy barrier for the segment L is U — max Ui . It 
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is important to stress here that we are not trying to find 
the optimal hop for the segment L and the associated 
optimal energy barrier, but we are simply discussing the 
distribution of all possible hops of the segment L from 
the initial to all available final metastable states. In order 
to construct all possible configurations of the advancing 
segment, we now imagine redefining the network of ele- 
mentary hops by choosing a new unit L\ > L c so that 
the section L is now composed of, say, n\ sub-blocks of 
size L\ and fluctuating energy barriers Uu. The set of 
all possible configurations of the excitation L as it vis- 
its all available metastable states forms a new cluster of 
elementary excitations or subclusters L\. Notice that 
the subclusters of length L\ may in general be multiply 
connected, but the argument given above will still apply 
upon redefinition of the elementary unit. The global bar- 
rier for the motion of the new cluster will again be the 
maximum of the energy barriers of each subcluster. The 
linear structure of the string ensures that the cluster of 
all possible configurations of L remains singly-connected 
at all levels of rescaling. It is precisely this linear topol- 
ogy that enables us to carry out this rescaling procedure. 
Repeated application of this procedure will allow us to 
list all possible configurations of the segment in questions 
and yield the limiting distribution that we seek. Such a 
distribution, if it exists, must be stable under the max 
operation and therefore belongs to the class of so-called 
extreme distributions [ fl5| . To briefly summarize the the- 
ory of extreme distributions, we consider a set of identi- 
cally distributed independent random variables Xj, with 
1 < i < n - in the present problem these are all the energy 
barriers of the elementary spins and/or block of spins 
composing the segment L. Let M n = max{Xi, X n }. 
If F(x) is the probability of the event Xi < x and V n {x) 
the probability that M n < x, it can be shown that in 
the limit of large n the probability distribution of the 
maxima can be approximated by an asymptotic form 
V(x), given by the solution of the functional equation 
limn^oo V n {x + a n ) = lim^oo F n (x + a n ) = V(x), where 
V(x) is defined for — oo < x < oo. The extreme distribu- 
tion we seek has the form V{x) — exp[— exp(— x)\ [fT5| . 
The functional form of the distribution must remain in- 
variant when the elementary unit of the cluster of all 
possible hops is redefined according to the procedure de- 
scribed earlier. Since ["P(x)] n = exp[— nexp(— x)], this 
requires a n — In n. It thus follows that for the problem 
of interest here the appropriate variable x is E — lnZ, 
where we measure E in units of U c and L in units of 
L c . The probability distribution of energy barriers for a 
given segment L is then Vl(E) — exp[— L exp(— E)]. The 
corresponding distribution density pl(E) = dVL{E)/dE 
is 

PL (E) = Le- E exp[-Lex P (-E)}. (4) 

The typical barrier of a segment of length L scales then 
as E ~ U c \nL. The global distribution density W(E) 



of energy barriers for the string is obtained by integrat- 
ing pl(E) over L with the proper weight ul describing 
the density of the segments of length L. At the critical 
point F = 0, the cluster-sizes distribution is nj, ~ L~ v . 
At finite forces F the same distribution will describe the 
cluster network up to L = L opt . For the usual percolation 
clusters v = 1 + d/df, with d the dimensionality of the 
medium and D < df < d the fractal dimensionality of 
the clusters network O] . For anisotropic directed perco- 
lation describing the depinning transition the exponent 
may be different and will be discussed elsewhere, but the 
relation v > 2 holds. By evaluating the integral we find 

W(E) ~e- { »- 1)E/u <. (5) 

The corresponding probability density of finding a wait- 
ing time r = r exp(_E/T), with r Q a microscopic time 
scale, is given by 

*(T)~T(T /r) 1+a , a={v-l)T/U c . (6) 

The distribution of energy barriers given in Eq. H must be 
cut off at the energy U (F) corresponding to the optimal 
segment L opt (F) since the segments of larger scales are, 
on average, sliding freely and do not participate in the 
activated dynamics. The notion of "waiting time" makes 
no sense for such excitations. The average waiting time 
controlling the motion is therefore 

/T max 
dr*(r)r ~ exp[(l - a)C/(F)/T], (7) 

where T max — t exp\U{F)/T]. Since the motion is con- 
trolled by the largest barrier, one arrives at the average 
velocity given by v ~ u opt /(r) ~ exp[-(l - a)U(F)/T], 
recovering the result of the scaling theory || . 

We now discuss some observable consequences of the 
result given in Eq. [| The vortex lattice in the mixed 
state of high temperature superconductors is probably 
the best system for the experimental study of differ- 
ent aspects of glassy dynamics, since by tuning param- 
eters such as the applied field one can probe different 
regimes. For a driven vortex lattice the avalanche-like 
low temperature dynamics just described will manifest 
itself in the spectral density of the density-density corre- 
lation function, 5(k = 0,t) = (Sn v (r, t)8n v (r, 0)}, which 
is accessible by Hall-probe measurements. The density- 
density correlation function can be expressed in terms of 
the correlation function of the displacement field u(r,t) 
as S^k, t) = k 2 (u(k, t)u(— k, 0)). On large distances 
displacements are uncorrelated in a glassy system and 
(w(k, t)u{— k, 0)) oc (l/fc 2 )/ u (k, t). The corresponding- 
spectral density J u {lo) can be expressed in terms of the 
distribution function of energy barriers W(E), 
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Using W(E) from @ we obtain a power law density noise 
spectrum, 



We see that at low temperatures, T < U c , the exponent 
a«l, and the noise spectrum of vortex density fluctu- 
ations is nearly the 1// spectrum. 

We now turn to discuss the high temperature region, 
where thermal fluctuations renormalize the pinning en- 
ergy barriers. The characteristic temperature Td p sepa- 
rating the high temperature and low temperature regions 
is defined by the self-consistent equation Td p — U c (Td p ) 
H. For elastic manifolds with D > 2, U C (T) grows much 
faster than T |Q| for T > T dp , and a < 1. The results 
obtained earlier are therefore still relevant, with suitably 
renormalized parameters. 

The situation is more subtle for the case of an elastic 
string, corresponding to D = 1. Above the string depin- 
ning temperature, TJ p , the minimal pinning energy of the 
string becomes of order T, i.e., U C {T) ~T |§. The def- 
inition of the depinning temperature used here follows 
from a scaling theory based on the postulate that the 
creep barriers scale in the same way as the fluctuations 
in the free energy, i.e., the string statistical mechanics is 
controlled by a unique energy scale Such a scal- 

ing theory leaves numerical constants undetermined. To 
preserve the idea that a single energy scale controls the 
dynamics we define T^ p as the temperature where first 
a(T) = 1. This definition avoids the introduction of 
a new unphysical energy scale depending on the fractal 
dimensionality of the cluster of the pinned states. For 
T > TJp the motion of the string on scales L < L opt is 
still governed by the waiting times distribution function 
of Eq. | with a = 1, i.e., *(r) ~ (1/t 2 ). The average 
waiting time in this high temperature regime is then 

f=T ln(T max /T ) ~F-». (10) 

The characteristic energy and length of the segment hop- 
ping with this waiting time are E ~ Tln(f/T ) ~ hiF 
and L = L c (t/t ), respectively. The typical transverse 
displacement u of a length L of string is determined by 
Cu 2 /L ~ FuL. Defining the string velocity as v — u/f 
and recalling that the activated velocity contains a pref- 
actor ~ F ensuring that flux flow is recovered for F > F c , 
we obtain 

v~(F/F c ) 2 -», T>Ti p . (11) 

In the high temperature region, T > T^ pJ the string ve- 
locity vanishes as a power law. This represents a kind of 
marginal glassy dynamics which may be relevant to the 
high temperature, low magnetic field region of type II 
superconductors. We conjecture that this marginal dy- 
namics corresponds to a new glassy state in the regime of 



single vortex pinning which can be referred to as critical 
glass. 

In conclusion, we have developed a microscopic model 
for the glassy dynamics of elastic manifolds. In our model 
the low-temperature creep is governed by a power law 
distribution of waiting times. The dynamics at low tem- 
perature, T < Td p , and long times is dominated by the 
optimal barrier U(F), corresponding to the maximally 
pinned configuration. As a result, the mean velocity of 
the manifold vanishes exponentially with vanishing driv- 
ing force according to v ~ exp[— U(F)/T], as obtained 
earlier by a scaling approach |jj . This form for the macro- 
scopic response arises naturally as a result of the cutoff 
of the algebraic time distribution at the maximal wait- 
ing time. We also conjecture a qualitatively different 
marginal glassy response, with v oc [(F/F C (T)] 1- ' 1 , for 
a single elastic string and T > T^ p . Finally, we discuss 
how the power-law distribution of waiting times controls 
the noise spectrum of the vortex density and can there- 
fore be probed experimentally 
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